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In this work we investigate a numerical model of a coaxial RF-driven plasma jet operated at 
atmospheric pressure. Due to the cylindrical symmetry an adequate 2-D representation of 
the otherwise 3-dimensional structure is used. A helium-oxygen chemistry reaction scheme 
is applied. We study the effect of secondary electrons emitted at the inner electrode as 
well as the inserted dielectric tube and discuss their impact on the discharge behavior. We 
conclude that a proper choice of materials can improve the desired mode of operation of 
such plasma jets in terms of materials and surface processing. 



1. Introduction 

Microplasmas operated at atmospheric pressure 
have recently gained high attention |1 1. A par- 
ticularly popular class of microplasmas are non- 
equilibrium radio-frequency driven plasma jets, 
originally proposed by Selwyn and co-workers in 
1998 1 2 1. Many groups have investigated their 
prospect for surface modification; the studied 
processes include etching of tungsten, deposi- 
tion and etching of silicon oxide, and cleaning 
of thermolabile surfaces from contaminants. In 
order to realize a jet concept for coating ap- 
plications, e.g. homogeneous thin film depo- 
sition from C2H2, von Keudell and co-workers 
developed a micro-scaled axis-symmetric radio- 
frequency driven plasma jet operated at atmo- 
spheric pressure (3J. The schematics of their de- 
vice is depicted in Fig. [l| A stainless steel capil- 
lary is inserted into a ceramic tube leaving an an- 
nular space of 0.25 mm between the outer radius 
of the capillary and the inner radius of the ceramic 
tube. The ceramic tube is finally surrounded by 
an aluminum tube. This outer tube is now driven 
through a matching network by an RF voltage at 
13.56 MHz. In contrast the inner capillary acts 
as the grounded counter electrode. A feed gas 
(helium or argon) is guided in the capillary (red 
arrow) and the annular space between the capil- 
lary and the ceramic tube (blue arrows), while a 
small amount of a reactive gas is additionally in- 
jected through the capillary. Besides proving the 
successful operation of this device as depositing 
tool the coaxial jet is investigated experimentally 
by means of phase resolved optical emission spec- 
troscopy (PROES). Four different modes of the 
discharge are distinguished in (4j. 
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Fig. 1: Schematic sketch and simulation domain of the 
coaxial jet. The inner capillary (white) is electrically 
grounded while the the outer aluminum tube (black) 
is driven by an RF voltage. The ceramic tube (gray) 
consists of aluminum oxide. 

In this work we study the coaxial jet by means 
of a two-dimensional numerical simulation. We 
take advantage of the axis-symmetric geometry 
and resolve the radial and axial directions of 
the jet. In particular, we focus on the role of 
the secondary electrons emitted at the ceramic 
tube and the grounded capillary in order to 
characterize the 7-modes of the discharge. 

2. Simulation 

2.1. Description of the numerical code 

The simulations of this study are carried out 
with the fluid dynamics code nonPDPSIM de- 
signed and realized by Kushner and co-workers. 
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For a detailed description of the code see publi- 
cation 1 5 and the citations therein. Here, we just 
briefly discuss the implemented equations and the 
underlying physics. 

The code nonPDPSIM simulates the dynam- 
ics of weakly ionized plasmas in the regime of 
medium to high pressure. It takes into account 
the physics and chemistry of charged particles - 
electrons with mass m e and charge — e, ions with 
mass rrij and charge qj - and of the excited as well 
as the ground state neutrals (mass mj). For all 
species j, the continuity equations (particle bal- 
ances) are simultaneously solved, where Tj is the 
particle flux density and Sj is the source and loss 
term, respectively: 



diy 
dt 



(1) 



The fluxes are calculated from the momentum 
balances in the drift-diffusion approximation eval- 
uated in the local center-of-mass system. Dj and 
Hj are the diffusion constant and the mobility (if 
applicable) of species j. Further, E is the elec- 
trical field, and v is the mass-averaged advective 
velocity of the medium: 



% = njv-DjVnj + ^ m E. 



(2) 



For the electron fluid, additionally an energy bal- 
ance equation is solved taking into account Ohmic 
heating and the energy losses due to elastic and 
inelastic interaction with the neutrals and ions as 
well as heat conduction, 



d /3 



dt \2 



n e T e = j-E . . . 
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(3) 



To capture the non-Maxwellian behavior of the 
electrons, all electronic transport coefficients (the 
mobility /i e , the diffusion constant D e , and the 
thermal conductivity K/q ) as well as the electronic 
rate coefficients in eqs. 

@ and (|]) are calcu- 
lated by the local mean energy method: A zero- 
dimensional Boltzmann equation for the electron 
energy distribution /(e) and the transport and 
rate coefficients is solved for the locally applicable 
gas composition and various values of the electri- 
cal field. The tabulated data are then consulted 
in dependence of the fluid dynamically calculated 
electron temperature T e . 



The plasma equations are coupled to a modified 
version of the compressible Navier- Stokes equa- 
tions which are solved for the gas density p, the 
mean velocity v, and the gas temperature T. The 
contributions to the energy equation from Joule 
heating include only ion contributions; the heat 
transfer from the electrons is included as a colli- 
sional change in the enthalpy. The scalar pressure 
p is given by the ideal gas law. 

Finally, the potential <3> is calculated from Pois- 
son's equation. (The code works in the electro- 
static approximation such that E = — V$.) The 
charge density on its right hand side stems from 
the charged particles in the plasma domain and 
from the bound charges p s at the surfaces. The 
coefficient e — e$e t represents the permittivity of 
the medium: 



-V • sV§ 



(4) 



The surface charges are governed by a separate 
balance equation, where a is the conductivity of 
the solid materials and the subscript s indicates 
evaluation on the surface: 



dt 



^^(-V.f i + ^-)-V.(a(-V$)) 



(5) 



The dynamical equations are complemented 
by an appropriate set of boundary conditions. 
Electrically, the walls are either powered or 
grounded. With respect to the particle flow, they 
are either solid, or represent inlets or outlets: 
The flow is specified to a given flux, while the 
outlet flow is adjusted to maintain the pressure. 
Finally, it is worth mentioning that the actual 
implementation of the equations poses some 
difficulties due to the vast differences in the time 
scales of the dynamics of the plasma and the 
neutrals. These difficulties are overcome by the 
methods of time-slicing and subcycling. 

2.2. Details of the simulation case 

The described model is employed to simulate 
the coaxial jet described before. Fig. [T] provides 
(besides the schematic sketch of the jet) the sim- 
ulation domain. We scaled down the length of 
the original experimental jet configuration in or- 
der to reduce the number of grid nodes but pre- 
served the original radial dimensions. In our case, 
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the capillary is 4.25 mm long, the grounded elec- 
trode is 5.25 mm and the ceramic tube is 6.25 mm 
long, respectively. The capillary has an inner and 
outer radius of 0.1 mm and 0.25 mm, the inner 
radius of the ceramic tube is 0.5 mm. The sim- 
ulation resolves the two cylindrical coordinates R 
and z. Due to the azimuthal symmetry of the 
simulation "three-dimensional physics", e.g. spot 
attachment of the discharge at a surface, can not 
be resolved. 

The device is operated in a pure helium atmo- 
sphere at 10 5 Pa with helium as feed gas and oxy- 
gen as the reactive gas. The feed gas in the annu- 
lar space between the capillary and the ceramic 
tube is injected with a flow rate of 3 slm, corre- 
sponding to a maximum advective velocity of 80 
m/s. To avoid turbulences of the gas flow at the 
tip of the capillary the flow rate of the gas in- 
jected in the capillary is set to maintain a similar 
advective velocity: Helium with an admixture of 
1 % oxygen as reactive gas is guided in the cap- 
illary with a flow rate of 170 seem. The outlet is 
controlled to maintain a constant pressure. 

We apply a helium-oxygen chemistry reaction 
scheme which is described in detail in |6|. 

To study the role and influence of the sec- 
ondary electrons on both surfaces - the grounded 
capillary and the ceramic tube - we vary the 
secondary electron emission (SEE) coefficient 7. 
Additionally, we weight 7 for each ion species to 
account for the different probabilities of emission 
of secondary electrons depending on the ion 
species hitting the surface. 

3. Simulation results 

We ignite the discharge numerically by seed elec- 
trons in the annular space between the capillary 
and the ceramic tube. The main ionization pro- 
cess during the numerical convergence changes 
from volume to secondary electron dominated ion- 
ization. This characteristic is found in the cor- 
responding experiment by increasing the applied 
voltage. In the converged state the discharge runs 
in the 7-mode as observed by the experiment. 

Following, two distinct discharge regimes are 
discussed with simulation results. First, the SEE 
coefficient for the ceramic tube is set to 7 cer = 0, 
while in the second setting it is set to 7 cer = 0.1. 
The SEE coefficient for the capillary is 7 cap = 0.1 
in both cases. Apart from that all other simula- 
tion parameters are kept constant for both cases. 
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Fig. 2: Electron density (cm -3 ) and ionization coef- 
ficient due to secondary electrons (cm -3 s -1 ) without 
secondary electron emission at the ceramic tube. 
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Fig. 3: Electron density (cm -3 ) and ionization coeffi- 
cient due to secondary electrons (cm -3 s -1 ) with sec- 
ondary electrons at the ceramic tube. 

In order to investigate the influence of the sec- 
ondary electrons emitted from the ceramic tube 
on the discharge, two aspects are characterized. 
As presented first in Figs. [2]and[3]the spatial pro- 
file of the electron density as well as the ionization 
due to secondary electrons is analyzed. This com- 
parison is then followed by a comparison (com- 
pare Figs, [i] and [5]) of the dominating chemically 
active species, namely atomic oxygen and ozone. 

In both scenarios the influence of the SEE at 
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Fig. 4: Distribution of reactive oxygen species (cm ) 
in the effluent of the coaxial jet without secondary 
electron emission at the ceramic tube. 
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Fig. 5: Distribution of reactive oxygen species (cm -3 ) 
in the effluent of the coaxial jet with secondary elec- 
trons at the ceramic tube. 

the ceramic tube is clearly evident. E.g., when 
comparing the electron density distributions a 
shift towards the ceramic tube is clearly visible 
with SEE at the ceramic tube. This strongly indi- 
cates that these emitted electrons do play a signif- 
icant role in the spatial distribution of the plasma 
discharge, bending it towards the surrounding ce- 
ramic tube. 

With respect to reactive species the secondary 
electrons emitted from the ceramic tube again 
play an important role. As can be seen in the 
spatial profile of ozone (compare Fig. [5]), which 
expands significantly further into the effluent of 



the plasma jet. 

These results suggest that the SEE from 
materials adjacent to the plasma in a coaxial 
setup has a significant influence on the spatial 
profile of the plasma discharge. In particular a 
proper choice of these materials can alter the 
discharge behavior with respect to the reactive 
species which finally determine the applicability 
in terms of materials and surface processing. 

4. Conclusions 

In this work we investigate a numerical model of 
a plasma jet in coaxial configuration operated at 
atmospheric pressure. We show that the emission 
of secondary electrons at the ceramic tube 
adjacent to the plasma discharge has substantial 
impact on the spatial distribution of the plasma 
as well as the reactive species, e.g. ozone. Further 
we conclude that a proper choice of materials can 
improve the desired mode of operation of such 
plasma jets in terms of materials and surface 
processing. 
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